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Abstract 



Recently Pluhar and Weidenmiiller [Ann. Phys. (N.Y.) 297, 344 
(2002)] showed that the eigenvectors of the matrix of second moments 
of embedded Gaussian unitary ensemble of random matrices generated 
by /c-body interactions (EGUE(/c)) for m fermions in N single particle 
states are SU (N) Wigner coefficients and derived also an expression 
for the eigenvalues. Going beyond this work, we will show that the 
eigenvalues of this matrix are square of a SU{N) Racah coefficient and 
thus the matrix of second moments of EGUE(A;) is solved completely by 
SU (N) Wigner-Racah algebra. 
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I. INTRODUCTION 



Interacting finite quantum systems such as nuclei, atoms, quantum dots, 
nanometer-scale metallic grains etc. are governed by hamiltonians of low (compared 
to the number of particles in the system) particle rank. Therefore, for many purposes, 
the random matrix models appropriate for these systems are embedded random ma- 
trix ensembles of /c-body interactions originally introduced, via nuclear shell model 
calculations, by French and Wong and Bohigas and Flores [1]. For a system of m spin- 
less fermions in N single particle states (we will use fermions throughout this paper 
and turn briefly to bosons at the end) the embedded Gaussian unitary ensemble of 
random matrices of A;-body interactions [EGUE(/i;)] is generated by defining the hamil- 
tonian H, which is given to be /c-body, to be GUE in the /c-particle spaces and then 
propagating it to the m particle spaces by using the geometry (direct product struc- 
ture) of the m-particle spaces. Just as EGUE(A;), the EGOE(A:) and other embedded 
ensembles are defined [2] . With m particle space dimension given by — (^^^ , one 
has the unitary groups SU{N), U{Nk) and U{Nm) with EGUE(/c) invariant under 
U{Nk) and the embedding in m-particle spaces is defined by SU{N); note that a 
GUE in m particle spaces is invariant under U{Nm) but not the EGUE(A;), k < m. 
Very early, using the so called binary correlation approximation, Mon and French [3] 
and later French et al [4] derived some analytical properties of embedded ensembles 
valid in the dilute limit (defined by (A^, m,k) —>■ oo, m/N —>■ and k/m —>■ 0). How- 
ever only recently rigorous analytical results, valid for any {N, m, k) are derived for 
these ensembles by Benet et al. [5] and very soon Pluhar and Weidenmiiller (here- 
after called PW) demonstrated [6] that these results indeed follow from considerations 
based on the SU{N) embedding algebra. With all the m-particle matrix elements 
being linear combinations of the /c-particle matrix elements (see Eq. (9) ahead), the 



2 



joint distribution for the matrix elements will be a multivariate Gaussian. Thus all 
the information about EGUE(A;) is in the covariance matrix or the matrix of second 
moments (Eq. (10) ahead). PW have shown that the eigenvectors of this matrix are 
SU{N) Wigner (or Clebsch- Gordon (CG)) coefficients and derived the expression for 
their eigenvalues using a duality relation for EGUE(/c). The purpose of this paper is 
to show that the eigenvalues can be written as SU{N) Racah coefficients and thus 
the matrix of second moments is solved completely by SU (iV) Wigner- Racah algebra. 
To this end results for SU{N) Racah coefficients given in [7] and [8,9] are used. We 
will start with some basic results given in PW paper. 



II. BASIC DEFINITIONS AND RESULTS 



Let us begin with m particles in N single particle states (unfortunately in PW 
£ is used in place of N but to keep the notations same as in our earlier papers 
[2,10], we use N). The single particle (sp) creation operator al for any i-th sp state 
transforms as the irreducible representation {1} of U{N) and similarly a product of r 
creation operators transform, as we have fermions, as the irrep {F} in Young tableaux 
notation. Let us add that a U{N) irrep {Ai, A2, • • • , Ajv} defines the corresponding 
SU{N) irrep as {Ai — AAr,A2 — Aat, • • • , Aat-i — Xn} with N — 1 rows (there are 
also other equivalent ways of defining SU{N) irreps given a U{N) irrep [6]). This 
U {N) SU (iV) correspondence is used throughout and therefore we use U (N) and 
SU (N) interchangeably. In PW, {I''} is denoted by fr and we will follow this notation 
from now on. With Vr denoting irreps (and other multiplicity labels) of the groups in 
a subgroup chain of U {N) that supply the labels needed for a complete specification 
of any m-particle state (for the purpose of the present paper the subgroup chain need 

r 

not be specified), the operator JJ a| and a normalized r-particle creation operator 

i=l 
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-^HfrVr) behave as the SU{N) tensors T'^''"'" and -^T^'-'"'' respectively. Using the 
composition formula, 

rpfmVm ^ ^ Cj^^^")^^^ 7'^'=^'= T"^^"^ , S ^ TTl - k (l) 

where C'/j^^'/^„^ is a SU{N) CG coefficient, a m-particle state \fmVm) — ^\fmVm) |0) 
can be written as a product of k and s — m — k particle states as, 

\fmVm) = (D"^ E ^\fkVk) \fsVs) C^f^Tu^^ (2) 

VkVs 

Some properties of the CG coefficients, used later in simplifications, are 

f-lfab^ab — ( _l\'l>{fa,fbJab) (^fab'"ab f-yfab'"ab _ r~lfabV^ 

^faVa fbVb V ^) ^fb^b faVa ' ^ faVa fb^b faVZ fbV^ 



(3) 



nfabVab ^ ( _'i)<P{fa,fb,fab) ^(Afc) fyfaVg _ 

^faVafbVb V ) ^ ^^^^ fabVabfbVE 

QfaVa _ qOO _ ( (jfab'"ab_ V — (j}ab'"ab_ 

00 faVa ' f^Va foj^ \J d[f ) ' /a"^/ /a^b /a'^S' 

E((jfab'Vab_ \* (jf'ab'"'ab_ _ ^ ^ ^ ^ 

Va ,Vb 

In (3), is a function that defines the phase for the 1 <-> 2 interchange in the CG 
coefficients, d{f) is the dimension of the irrep / and / is the irrep conjugate to /. 
For fr = {1^}, fr — {l^"*"} and it also contains a phase factor as given in Eq. (53) of 
[7] (this is also seen easily by comparing the second and third equalities in (3) with 
the corresponding relations for the standard CG coefficients for angular momentum 
[11]). Similarly with fa = fk, fb = fk, one has 

fab^g. = {2''l''-"'}, u^O,l,...,k. 

Note that go = {0} for SU {N) and also g^, = y^. The function 4> in Eq. (3) is of the 
form 0(Ai, A2, A3) = F{\i) + F(A2) + -^(Aa) where F is some function with F(A) an 
integer and F{\) = F{\); for SU{N) irreps that appear in this paper these results 



are valid as can be seen from Eq. (60) of [7]) except that there can be an overall N 
dependent factor which will not change any of the final results. 

SU{N) irreducible tensors Bk{gu^v) constructed out of A'^{fkVk)A{fkv'f?) are de- 
fined by, 

Bu{9.uj.) = E A\hv,)A{hv',) Cl^j^^ (4) 

It is useful to note that the tensors S's in (4) multiplied by k\ are, to within a 
phase factor, same as the tensors defined in Eq. (48) of [7] . Wigner-Eckart theorem 
decomposes the matrix elements of Bk{g^uj^) in m-particle spaces into a reduced 
matrix element ( 1 1 1 1 ) and a CG coefficient, 

{fmVm I Bkig^UJ^) I fmV'm) = (/m 1 1 Bk{g^) \ \ fm) C^.^""" -p-^ (5) 

Two important properties of Bk{g^u!,^) are, 

{{Bk{g^uj^))f = //Vfc5p^,{o} 

(6) 

{{Bk{guu;^)Bk{g^Lj^)))'' = 6g^g^ 6^^^^ 

In (6), (O)'^ denotes trace over the /c-particle spaces. The first equality in Eq. (6) 
easily follows from the fact that here g^ = {0} as traces are scalars with respect to 
SU{N) and then applying the fifth equality in Eq. (3). Similarly the second equality 
in Eq. (6) follows from the fact that only g^, — g^ will give a scalar. Now we will turn 
to EGOE(A;). 

Consider a /c-body hamiltonian 

H{k) = E V,av>,{k)A\f,Va)A{f,v,) (7) 

where Vv^y^^{k) are matrix elements of H{k) in /c-particle space and form a GUE, i.e. 
^vaVbi^) are independent Gaussian variables with zero center and variance given by 
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In Eq. (8) the overline indicates ensemble average and is ensemble averaged vari- 
ance of H{k) in /c-particle space. The m-particle matrix elements oi H{k) are, with 
s — m — 

H.wjk) = {fmvi\H{k)\f^vi) = {-) y: (^/r;iO* (9) 

Va,Vb,Vs 

Eq. (9) is obtained easily by substituting the definition (7) for H{k), then inserting 
complete set of states between and A operators and applying Eq. (2). The 
EGUE(A:) in m > k spaces is defined by Eqs. (7)-(9). Now it is clear that for 
any analysis of EGUE(A;) all one need to know is the covariance between any two 
m-particle matrix elements Hyi (k) and this defines the matrix of second moments, 



Ayl ^4 .^3 y2 — Hyl y1 (k)Hy3 ^4 ( /C ) (lO) 

As stressed by PW, most important step in EGUE(A;) analysis is to derive a "gen- 
eralized eigenvalue expansion" of A defined by Aij = J2k C-ikEkCjk with the 
eigenvalues and Cjk the eigenvectors such that Ek are positive and C's hermitian. To 
this end, it is useful to consider the unitary decomposition of H{k) in terms of the 
SU{N) tensors Bkig^tOu), 



H{k) = J2 Bu{g.uj,) Wg^^Sk) (11) 

The expansion coefficients Wg^^^{k) are easily given by 

Wg^u^Ak) =« H{k)Bk{9.uj,) (12) 

and this follows by using the definition (11) and Eq. (6). Most significant property of 
the W coefficients is that they are independent Gaussian variables with zero center 
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and variance given by (derived using Eqs. (11), (4), (7) and (8) in that order and 
using the orthonormal properties of the CG coefficients) , 



W,^.Sm,,.,{k) = — S,^,^ 5^^^^ (13) 



III. MATRIX OF SECOND MOMENTS 



Firstly we will derive an expression for the covariance Hyi^^2^{k)H^-A^y4^{k) in terms 
of SU {N) CG coefficients and Racah coefficients and then turn to the eigenvalues and 
eigenvectors of A, the matrix of second moments. Applying Eqs. (10), (11) and (5) 
in that order gives. 



Hyi ^2 ik)Hy3 (k) 



E (/rn^m I Bkig^^u) Wg^u^Ak) I fmVl) {fmVl \ Bk{g^UJ^,) Wg^^^k) I fmVt 



TT E \fmVln I BkiQvU^u) I fmVl,) {fmvl, \ Bk{g^UJ^) \ fmV, 

gv^^v 

y2 2 a; 

^ E \{fm\\ Bk{gu) W fm)\ ^-rr 3 

1^ gvijJv 

(14) 

Eqs. (55) and (56) of [7] together with Eqs. (3) (see also the remark just after Eq. 
(4)) allows one to write the reduced matrix element in Eq. (14) as a SU{N) Racah 
or [/-coefficient, 

2 (Nm)'^ 2 

|(/m II Bk{gu) II fm)\ = -Ji TTTT 7 \U{fmfN-kfmfk]fm-kgv)\ (15) 

(^\gu)V^m-k) 

With {2'^l^-2^} = {1^} (8) {l^-'^} - {l'^-^} ® {1^-^+1} where ® denotes Kronecker 
product, the dimension d{gy) is given by 
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Before going further let us define, for a given (m, N), a function A''(/c), 

m — f \ / N — m -\- k — v 



A'^ik)^^--) : j (17) 

Now, using Eq. (16) and substituting the expression given by Eq. (61) of [7] for the 
[/-coefficient in (15), it is seen that 

|(/m||5fc(^.)||/m)|' = A"(m-A;), i/ = 0, 1, 2, . . . , A; (18) 

Combining (14) with (18) yield an expression for the covariance between H matrix 
elements in m-particle spaces. 



Hyi y2 (k)H^3 y4 ik) 



_ ^ {N^f{^r ,rrf, . .... 

E - k)} c;-; c;-; 

iVl. „ , , JmVmJm.Vm JmlJjnJmUm 

l'=U,l,...,K;U)if 

(19) 

Eq. (19) will be useful in deriving expressions for the moments of H spectrum, i.e. 



(Hp)"^. However the disadvantage of (19) is that it is not in a proper form to give 
the eigenvalues and eigenvectors of the matrix A in Eq. (10). In order obtain them, 
the CG coefficients in Eq. (14) should be changed to Cf^^—^ Cf'^>; and this 
can be accomplished by a SU{N) Racah transform. Using Eq. (3.2.17) of [9] one has 
for example. 

Then Eqs. (20) and (3) will give. 



8 



fm "^m fm '^m fm "^m fm '^m 



J/ \ ^ yjmjmjmjmiyvyn)'^ f ,,i f ,a ^ f ,,3 f ,,2 



(21) 



Finally Eqs. (21) and (15) combined with (14) produce the generalized eigenvalue 
expansion of the matrix of second moments A, 



A. 



1) IJ ■ '1)^ 



H,i ,2 ik)H,, ,4 ik) = y cf 4V ^ ?; ^ 



E 



\U {fmfN-kfmfk'i fm-k9u)] U {fmfmfmfm'i 9v9ti) 



(22) 



Obviously the quantity in the curly brackets in Eq. (22) gives the eigenvalues of A 
and the C's are eigenvectors. 

IV. EIGENVALUES AS SU{N) RACAH COEFFICIENTS AND THEIR 

APPLICATIONS 

In order to proceed further, it is useful to consider 6j symbols of SU (N) and they 
are defined by (see Eq. (3.2.18) of [9]), 



[/(AiAsAAa; A12A23) = V^(Ai2)d(A23)(-l)'^(^^'^^'°)+'^(^^^'^^'^)+'^(^^'^^'^^^) 



Ai A23 A 
A3 A12 A2 
(23) 



In (23) A's are SU{N) irreps and the four couplings involved in the [/-coefficient 
are assumed to be multiplicity free (for the applications in the present paper this 
assumption is always valid) . Symmetry properties of the 6j-symbol appearing on the 
r.h.s of (23) are well known [8,9]. In the present analysis, the Biedenharn-EUiott sum 
rule extended to SU{N) [8,9] plays a central role. This sum rule relates a product 
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of three Racah coefficients (weighted appropriately by dimension factors and phase 
factors with the irreps in the Racah coefficients appearing in some particular order) 
with sum over a common irrep label to a product of two Racah coefficients. After 
converting the Racah coefficients in (22) into 6j symbols of SU {N) using Eq. (23) 
and then applying the symmetry properties of the 6j symbols, it is seen that the 
sum in the square brackets in Eq. (22) is exactly in the required form. Applying the 
Biedenharn-EUiott sum rule, the sum then simplifies to 

U {fmfN—m+kfmfm—k > fk9ti) 



Nk dig,,) 

Here = 0,1, . . . ,m — k. Now the eigenvalues of the matrix A, in terms of the U (N) 
Racah coefficients is given by 

\2 /AT \2 /m\2 

= m VaA [U{fmfN-m+kfmfm-k; IkQ,)? ] = 0, 1, ■ ■ ■ , 771 - k (24) 

-'Vfe a[gij,)[l\k) 

with degeneracy d{gij,) (see Eq. (16)). Eq. (24) is the central result of this paper. 
With this, the matrix A is completely specified by the U{N) Wigner and Racah 
coefficients. Now substituting the formula (Eq. (61)) of [7]) for the [/-coefficients in 
Eq. (24) produces the result of PW, 



A^l ^4 .^3 ^2 = H^l ^2 {k)H^3 ^4 {k) = V Cf _ Cf „ E. (25) 

where, 

E^^^K^^ik), pi^O,l,---,m-k. (26) 

^^k 

Note that the function A'*(A;) is defined by Eq. (17). 

Information about EGUE(/c) is contained in the ensemble averaged moments Mp = 



(Hp)"^ and the bivariate moments T,pq — (Hp)"^ {H^)^. In deriving the formulas for 
the lower order moments, we will show the usefulness of Eq. (19). Obviously, ensemble 
averaged centroid is zero and the variance is 
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(27) 



The second equality follows from (25) and the final result follows by applying (3); 
note that C!'''^''^-^ = VN^^uo- The variance in (27) is in X^/Nj, units and 
this factor is dropped as all the quantities we consider from now on are all scaled 



with respect to {(i/^)'"}^/^. As the third moment is zero, we will turn to the fourth 
moment, 

1 



i V 7 



N ^ 



+ 



fm 



X 



E 



y {A"(m-A;)} C'^^, _— C 



'"^ „i „J ,,fe' „i I J/=0,l,...,A;:a;,y 



X 



^H=U,l,...,m—K;(jj^ 

-t min{k,m—k} 

= 2 [A0(^)]^ + — Y: K^{m - k) K^{k) d{v) 



(28) 



The second equality in Eq. (28) follows by applying Eqs. (11) and (13). In the third 
equality, it is easy to recognize the first term. The second term follows by applying 
Eq. (19) and (25) to the two pieces in the corresponding term in the second equality. 
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The final result follows by applying the orthonormality of the CG coefficients. Eqs. 
(27,28) will give the excess (72) parameter of the density of eigenvalues of EGUE(/c), 



72 



{my 



rmn{k^-k} _ ^^^^ 



v=0 



(29) 



Turning now to the lowest bivariate moment En, it is easily seen that 

1 ^ ^ ^ 1 



En = {Hr W 



Nrr 



A°(m - k) 



(30) 



Applying (19) and recognizing that only u — will contribute to the traces give 
immediately Eq. (30). However Eq. (25) generates a different formula for En and 
equating it to (30) gives the identity (derived in PW using the duality transformation), 
— — ^ A''{k)d{u) — A°(m — k). The variance of the distribution of centroids of the 
H spectra over the EGUE(A;) ensemble is 

En 1 AO(m - k) 



En = 



N„ 



AO (A;) 



(31) 



Finally E22 is given by 



E22 = {H^r{H^y 



E 





2 1 |2 


Ha j 


111 k' I \ 



H 



11 h' 



+ 



"^m f'^m ■I'^rn ^'^rn 



E 



(32) 



[K\k)f + 



m—k 



u=0 



Here in the second equality used is the property x'^y'^ = (x^) (y^) +2[xy]^ of Gaussian 
variables x and y. Similarly the final result follows by applying (25) to the second 
term in the second equahty and simplifying the CG coefficients using Eq. (3). Now, 
the variance of the distribution of the variances of the H spectra over the EGUE(/c) 
ensemble is, 

^ E22 - 2 



-'22 



- 1 



[A^(fc) 



AO(A;) 



(33) 
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Eqs. (27,29,31,33) are also given by Benet et al [5]; this paper neither gives details of 
the derivations nor uses SU{N) Racah coefficients. Also in this work. En, E22 and 
72 + 1 are denoted by S, R and Q respectively while we have followed Ref . [3] . 



V. CONCLUSIONS 

Going beyond PW, matrix elements of the matrix of second moments are written 
explicitly in terms of SU (N) Wigner and Racah coefficients and this result is obtained 
by recognizing that the reduced matrix elements of Bk{gv) are SU{N) Racah coeffi- 
cients. With this one has Eq. (19) and this is converted into the generahzed eigenvalue 
expansion form by first applying a SU{N) Racah transform and then applying the 
Biedenharn- Elliott sum rule extended to SU{N). This gives the eigenvalues of the 
matrix of second moments exphcitly in terms of SU{N) Racah coefficients (Eq. (24)). 
The two different forms given by Eqs. (19,25) for the covariances of m-particle H 
matrix elements, give in a simple manner the formulas for the low order moments Mp 
that define the state density and the bivariate moments Epg that give information 
about fluctuations. 

Although EGUE(A;) for only fermions is considered in this paper, all the results 
in fact translate to those of EGUE(A;) for bosons by using the well known N —N 
symmetry [12,13], i.e. in the fermion results replace N by — and then take the 
absolute value of the final result. For example, the m boson space dimension is 



d(m) = 



= {{N — m + l)m}- More importantly the eigenvalues of the matrix 



of the second moments are 



-N — m -\- k — ly 
k 



m — ly \ I N -\- m -\- ly — 1 
k I \ k 



(34) 



This result was exphcitly derived in [14]; see Eq. (14) of this paper. Moreover for 
bosons, {k} O {k^'^} = {2i/, i/^-^} , i/ = 0, 1, ■ ■ ■ , /c. Also, the N ^ -N 
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symmetry and Eq. (16) give d{giy) = {{N + u - 1)^}^ - {{N + i/ - 2),,_i}^ and this 
is same as Eq. (15) of [14]. Similarly Eqs. (27,29,31,33) extend directly to the boson 
EGUE(A;) with A'^{k) replaced by Ag(/c) defined in Eq. (34). In addition, for fermions 
to bosons there is also am ^ N symmetry and this connects fermion results (say for 
Mp and Ep^) in dilute limit to boson results in dense limit [13]. 

Recently there is considerable interest in mesoscopic physics to study EGUE(/c) 
for fermions with spin [15] and here the embedding algebra is U{2N) D U{N) ® 
SU{2) with SU{2) generating spin. The approach presented in Sections II-IV is 
being apphed to this system; some useful results for the U{2N) D U{N) ® SU{2) 
Wigner-Racah algebra are available in Refs. [16,17]. Finally, Wigner-Racah algebra 
analysis of embedded ensembles with more general group symmetries (see [2,6,17,18] 
for examples) should be possible in future, thus opening up a new direction in random 
matrix theory. 
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